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ID diagonally disordered chain with Frenkel exciton and long range ex- 
ponential intersite interaction is considered. It is shown that some states of 
this disordered system are delocalised contrary to the popular statement that 
all states in ID disordered system are localised. 

I. INTRODUCTION AND MAIN RESULT 

It is well known that all wave functions of translationary symmetric systems are delocal- 
ized. One of the most interesting properties of the homogeneous disordered systems is the 
possibility of localised wave functions. 

The mathematical problems of the theory of disordered systems are very complicated 
and for this reason the theory of disordered systems is not so well developed as the theory 
of symmetric systems. Despite this fact some statements related to disordered systems 
are considered to be well established and reliable. The above mentioned occurrence of the 
localised band is one of them. The next example of statement of this kind is that all states 
in ID disordered system are localised [5]. Recently appeared the reports [3,2] about the 
delocalisation in ID systems with intersite interaction in the form: J n ^ m = J/\n — m\ u , 1 < 
v < 3/2. In this letter we consider ID diagonally disordered chain with exponential intersite 
interaction and present arguments (computer simulations and theoretical treatment) in favor 
of partial delocalisation in this system. In this section we describe the system and present 
numerical results and in the next section we review the reasons which made us to study this 
system and present the approximate expression for the mobility edge. 

Let us consider ID Frenkel exciton in diagonally disordered chain. The mathematical 
problem is redused to the following random matrix of the Hamiltonian: 

H r y = e r 5 r y + w(r — r') r, r' — 1, N, (1) 

where 

w(r) — Vo exp — \r/R\, (2) 
Random values e r are supposed to be independent and having the distribution function: 

p(e) = | £ ^ ^ (Anderson's model) (3) 
I other cases 
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The thermodynamic limit JV — > oo is implied. To separate the localised and delocalised 
wave functions one should use some criterion of localisation. We use the number of sites 
covered by the wave function [8] determined as follows. Let us consider some eigen vector 
^ of the Hamiltonian (1) with components ^ r ,r = 1, JV. What contribution one should 
ascribe to the arbitrary site r? It is naturally to accept that this contribution is zero if 
\^ r \ 2 = and equal to unit if \^ r \ 2 =max{|\I/i| 2 , |\I/Ar| 2 }. So we come to the conclusion 
that the contribution of the arbitrary site r is |\I / r | 2 /| , I / | 2 na . r . The total number of sites n(\I / ) 
covered by normalised eigen function \& is the sum of contributions of all sites: 



n(¥) = £ 
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Delocalisation in (1),(2),(3) appear when R » 1. Below we study the properties of the 
eigen vectors of (1) with R = 20, v = 0.5, A = 4, JV = 1000. 

The dependance of number of sites covered by the wave function against corresponding 
energy for the Hamiltonian (1) is presented on fig.l (top). It is seen that n(E) is drastically 
increasing for energies higher than Eq ~ A. Additional calculations shows that n(E) do not 
depend on the number of sites JV for E < E and is ~ JV for E > E . For all these reasons 
we conclude that states below E are localised and states above E are delocalised. 

II. QUALITATIVE TREATMENT 

On our opinion the main properties of the above model which are responsible for the 
delocalisation are long range of intersite interaction R and the fact that function p{e) differs 
from zero only in the finite region. For these reasons for the qualitative interpretation we 
apply the following exactly solvable simple model of disordered system. Let the radius of 
interaction be infinite and write down the simplified Hamiltonian in the form: 

v 

H ry = 5 ry s r + — , r, r' — 1, JV (5) 

Taking advantage of the coherent potential approximation [6,5] one can show that the 
density of states for Hamiltonians (1) and (5) is coincide in the limit R — > oo, v — > 0, 2Rv = 
v. We show below that the Hamiltonian (5) has one delocalised and JV — 1 localised eigen 
functions. Consequently at least one delocalised function should appear in the set of eigen 
functions of the Hamiltonian (1) in the limit R — > oo. The desire to see how this take place 
was the starting point for our study of the Hamiltonian (1) with R » 1. Now let us turn 
to the proof of the above properties of the Hamiltonian (5). 

The equation for eigen vector e and eigen value A of the Hamiltonian (5) can be written 
in the form: 
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5 = Ee r (6) 



(6) gives an explicit expression for the eigen vectors of (5) as a functions of r and eigen 
number A. By substituting e r in the formula for S one can obtain the equation for the eigen 
values A: 

1^1 1 

^g— -f(A) = - (7, 

For Anderson's model (3) all quantities e r are differs from each other. For the graphical 
treatment of (7) the qualitative form of T(E) function is presented on fig.2. From fig. 2 one 
can see that N — 1 eigen values are belong to [0, A] and the last eigen value E m is not belong 
to this interval and in the thermodynamic limit can be determined from the equation: 



v ; J E-e A E-A v 

whence 



E m = z (9) 

1 — exp —(A/v) 

So in the thermodynamic limit E m is separated from any of e r by finite interval. From 
(6) one can see that sharp extremums of the wave function related to the localisation can 
appear if A G [0, A]. E m do not belong to this interval and we come to the conclusion that 
the corresponding eigen vector in the case of homogineous disorder is delocalised. Now let 
us show that all others eigen vectors are localised. For this reason introduce the Green's 
function in t-representation exp(?tH) r r which describe the dynamics of the wave function 
on the site r if it was equal to 1 on this site at t — 0. If the finite part of eigen states of the 
Hamiltonian H is delocalised this function goes down to zero when t — > oo and N — > oo. If 
Green's function do not decrease it means that the main part of eigen vectors is localised 
and the part of delocalised states is extremely small [5]. It is convenient to introduce the 
Green's function in E-representation: 

expitH= f G(E - i6) exp(iEt)dE 5^+0, t>0 (10) 

In the case of Hamiltonian (5) the Dyson's series for G: 

11111 11 

G rr , = b rr < + — v + — v 2 T(E) + (11) 

E — e r N E — e r E — e' r N 1 1 E - e r E - e' r 1 ; 

+i« 3 r(£) 2 - — + ... 

N 1 ' E-e r E-e' r 

can be exactly summed and give the following expression for G rr : 
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where 

1^1 1 V 

g(E) = i-y——^T(E)^±-\n-f— (13) 
yv y Nf^E-Ei A E — A y J 

In the thermodynamic limit the term ~ 1/N should be omitted and we come to the 
conclusion that the Green's function have a single pole E = e r . This corresponds to the 
oscillations of the wave function with constant amplitude and we can conclude that the main 
part of states are localised. It is easy to see that above described oscillations corresponds to 
the wave function localised on the site r and having an eigen value e r + 0(1/N). From fig. 2 
one can see that there are N — 1 eigen values of this kind and we come to the conclusion 
that the Hamiltonian (5) have N — 1 localised states and one delocalised with eigen number 
E m (9). 

Note that the appearance of the separated delocalised state for (5) is possible only if 
the distribution function p[e) is differ from zero in finite interval. For this reason we expect 
that delocalisation in (1) is also possible if p(e) is differ from zero in finite interval or at 
least goes down to zero rapidly enough. This statement confirms by calculations for Lloyd's 
model with p(e) = (l/7r)A/(A 2 + e 2 ) when no delocalisation was found. 

The energy dependance of number of covered sites for the Hamiltonian (5) is presented 
on fig. 1 (bottom) for v = 20. Other parameters are the same as for the top picture. One 
can see that finiteness of the interaction radius results in appearance of delocalised states in 
the gap [A, E m ] but the boundary energy of spectrum and the mobility edge are the same 
for both Hamiltonians (1) and (5) and are equal E m (9) and A respectively. 

III. EXPLANATION 

We explain the above numerical results of section 2 by studing the relative values of gaps 
between energy levels of undisturbed system and shifts produced by disorder [4]. Let us 
consider the spectrum of undisturbed system i.e. when H = W where matrix W has the 
following elements W rr i = t> exp(— |r — r'\/R). Its boundary energies can be calculated as: 

sh(lAR) 

u± = v - 



ch(l/ J R) ± 1 

If R » 1 this spectrum displays a very narrow peak at the low boundary energy u + so the 
system looks like degenerate. In this region these quasidegenerate states are strongly mixed 
by a disorder and the problem can be solved in a manner similar to that of previous section 
and the spectrum coinciding with p(e) can be obtained. Note that corresponding wave 
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functions are localised (see the previous section). Now let us turn to the dilute part of the 
spectrum. When R » 1 one can say that, in fact, the "dilute " part of spectrum occupy the 
interval from u + up to w_ excluding the above mentioned very narrow peak at E — u + . But 
this narrow peak in presence of disorder give rise to the interval of strong mixing [u + , u + + A] 
with localised states. For this reason the interval [u + ,u+ + A] must be excluded from the 
"dilute " part of spectrum. In the "dilute " region of spectrum twice degenerated levels of 
undisturbed system | ± q > are mixed by disorder 5 rr re r . The diagonalization of disordered 
part of the Hamiltonian 5 rr ie r in this two-dimensional space | ± q > r = (l/v^/V) exp(±igr) 
gives new eigen values: 

^± = fq + ^0 ± \£q\ 

where 
and 

fq = v J2ex.p(iqr - \r\/R) 

r 

- undisturbed spectrum of W. The direct calculation shows that the increment of undis- 
turbed spectrum e =t \£q\ is a random value with main value and dispersion ~ l/y/N. Here 
we neglect the constant shift A/2. For the gaps between levels of undisturbed system be 
greater than this increment (weak mixing) the following conditions must hold: 

1 A77 3 



Na(E) v^V ' Wl 
where cr(E) - is matrix W normalised density of states which can be calculated as: 

a(E) = 1 , V0SH1/R) V = votHl/R) 

71 EtJe 2 -(V - Efdi 2 {l/R) 

Energy Eq at which this density gets on its minimum can be calculated by formula: 



3 + ^9 - 8th 2 (l/i?) 
E ° = 4th(l/i?) 

At this energy the gap between undisturbed levels gets on its maximum and if it is 
greater than A^/v^/V one should expect weak mixing and extended (having size comparable 
with systems size) states. Equalising this two quantities allows one to estimate the system 
size iVo at which an extended states still take place: 

,a(E )r)A) 

When R » 1 this formula can be simplified: 

AW.687(^)V 



Nn 
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When N >> N one should expect strong mixing in the region of dilute spectrum (low 
density of states) and (may be and may be not) localisation. Anyway the size (4) of states 
in the region of E can be much greater than that at E ~ A. So in the region E G [u + , u + + A] 
systems states are originate from quasidegenerate part of matrix W spectrum, strong mixing 
take place and all states in this region are strongly localised. In residual part of spectrum 
when above conditions are satisfied mixing is very low and corresponding states can be 
extended. In the above numerical experiments R = 10 A = 4 and v = 1, so iVo ~ 10000. In 
this case systems with sizes iV ~ 1000 display the "mobility edge" at E — A where energy 
levels of disordered system originated from quasidegenerate part of undisturbed spectrum 
are adjoin to the region of low mixing with extended (under the above conditions) states. 
In our opinion even when N >> Nq the states in the region E G [u+ + A, u_] will be much 
"more extended" than those in region E G [u + ,u + + A] and some illusion of mobility edge 
at E = A still can appear. 
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FIG. 1. Energy dependance of number of sites covered by the wave function for the Hamiltonian 
(1) with N = 1000, A = 4, w(r) = v exp-\r/R\, v = 0.5 , R = 20 (top). The same for the 
Hamiltonian (5) with v = 2Rvq = 20 (bottom). 
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FIG. 2. 



